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Abstract 

We show that the particles in the Calogero-Sutherland Model obey fractional 
exclusion statistics as defined by Haldane. We construct anyon number den- 
sities and derive the energy distribution function. We show that the partition 
function factorizes in the form characteristic of an ideal gas. The virial ex- 
pansion is exactly computable and interestingly it is only the second virial 

coefficient that encodes the statistics information. 
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There has been some progress recently in understanding systems of particles with 
fractional exclusion statistics as defined by Haldane ||. Motivated by physical examples, 
such as quasi-particles in the fractional quantum Hall systems and spinons in antiferro- 
magnetic spin chains, Haldane formulated a generalized exclusion principle as follows: He 
considered many particle systems with finite dimensional Hilbert spaces where the dimen- 
sion of the single particle Hilbert space depends linearly on the total number of particles 
present. Then the exclusion statistics parameter g is defined by Ad = —gAN, where Ad is 
the change in the dimension of the single particle space and AN is the change in the number 
of particles. Thus g is a measure of the (partial) Pauli blocking in the system - g = 0(1) 
corresponds to bosons (fermions). In a recent work we have generalized the definition 
of g to the case of infinite dimensional Hilbert spaces so as to apply Haldane's definition to 
systems in the continuum with no sharp cutoff in energy. If for such a system a virial ex- 
pansion exists and all the virial coefficients are finite in the high temperature limit, then the 
exclusion statistics parameter g is completely determined by the second virial coefficient in 
this limit. In many systems it is also possible to relate g to the exchange statistics parameter 
a. 

In a remarkable paper [F2J, Wu has recently considered a system with flat dispersion and 
assuming that Haldane exclusion principle holds (as in the case of anyons in a magnetic 
field in two dimensions confined to the lowest Landau level) has shown that the statistical 
distribution function of this system is, 



Here e is the single particle energy, and w(e) is given as the solution of the equation, 



where (3 is the inverse temperature and \i is the chemical potential of the system. The 
distribution function smoothly interpolates between the Bose (g = 0) and Fermi(g = 1) 
distribution functions. It therefore seems likely that the distribution function discovered by 
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Wu may be valid in a more general context. Indeed it could be the general form for the 
ideal (exclusion) anyon Q distribution function. We will demonstrate that this is so for a 
model hamiltonian system in one space dimension with inverse square law potential- the 
Calogero-Sutherland Model(CSM) §. 

It has lately been recognized that models with the inverse square law interaction, the 
CSM H, Haldane-Shastry |§ model and other related models behave very much like ideal 
gases and that the particles have fractional statistics 0. In this paper, we concentrate 
on the CSM, which is a model defined in the continuum. Though the spectrum of the 
model and its thermodynamics are exactly known || the ideal anyon gas nature has not 
been brought out so far. To this end we apply our definition of exclusion statistics |L| to it 
and show that the particles have fractional exclusion statistics. Next we define the anyon 
number density in the energy space and prove that, in the thermodynamic limit, its average 
is given by the distribution of the form discovered by Wu. We then show that the partition 
function factorizes into factors obtained by integrating the distribution function with respect 
to the single particle energies. This clearly demonstrates the ideal gas nature of the system. 
Finally, we comment on an interesting feature of the exact equation of the state. 

We work in the fermionic basis of the CSM. The hamiltonian of the system of interacting 
fermions in this model is given by (H = 1), 
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where the particles are confined in a harmonic well and the thermodynamic limit is obtaining 
by taking uj — > 0. The particles can also be put on a circle with Xj\ replaced by the chord 
length \sin^(xi — Xj)\ in the interaction term. This model has the same thermodynamic 
limit as the one in eq.([|). At special values of the coupling the model can be mapped on 
to particular matrix models ||. It has recently been shown that a system of 2-d anyons 
(defined through the exchange phase) confined to the rim of a disc can also be mapped on 
to the model on a circle ||. 

The spectrum of CSM is exactly known. The states can be labelled by a set of fermionic 



occupation numbers {n k }, k = 1, oo, = 0, 1. The energy is given by, 

£[K}] = E - ^(i - 0)-^ — -> ( 4 ) 
fc=i 

where e k = ku, N = J2T=i n k- We note that this spectrum is identical to the spectrum of 
quasi-particle states in a gaussian theory of compact bosons [flj with radius R = l/^J\g) 
with the following identification- while k in the above equation is a state label in CSM, 
in the gaussian theory of compact bosons k refers to the box quantized momenta. The 
following discussions thus applies to this case also. As can be seen from eq.(|), the effect of 
the interaction is that each particle shifts the energy of every other particle by a constant 
uo(g — 1). The energy functional can also be written as, 

oo oo 

E[{n k }} = J2e k n k -uj(l-g) n fci n fe2- ( 5 ) 

fc=l kx<k 2 =l 

The exact spectrum of the model is thus reproduced by an effective hamiltonian of quasi- 



particle with constant density of states and constant Landau parameters flQ |. As we had 
discussed in our earlier paper this scale invariant energy shift is the basic reason for the 
occurrance of nontrivial exclusion statistics. We had also shown that a spectrum of the form 
in eq.(§]) results in the exclusion statistics parameter being equal to g . 

Thus the particles in the CSM are anyons in the sense of Haldane with the exclusion 
statistics parameter equal to g. We will now define anyon number densities and show that 
their thermal average is exactly given by the distribution function in eq.(|I]). To this end, we 
define N(e, 0) as the number of particles with energy e k < e, ie., N(e, 0) = J2tLi #( e — £k) n k, 
where 9{x) = for x < and 1 for x > 0. We may now define the shifted single particle 
energies as 

e A {k,g) = e k -u(l-g)N(e k ,0). (6) 

We identify this shifted energy with anyon single particle energy since the total energy can 
now be written as 

oo 

E [i n k}] = e ^(k, g)n k (7) 

k=l 



The number of particles with energy less than e is then N(e,g) = J2kLi — e A(k, g))^. 
The anyon number density in the thermodynamic limit is then given by 

„, (£ , 9)=limlim ^i±^|)^M, (8) 

where the limit uj — > is to be taken first. 

We will now derive a differential equation for N(e, g) and solve it to obtain the anyon 
distribution function. Consider a state with N-particles labelled by i = 1, N ordered 
such that ki + i > fcj. Thus for this state iV(e^,0) = i — 1. The i th shifted energy is 
£A(ki,g) = £k t -u(l-g)(i - 1). Note that e A (k i+1 ,g) - e A (k h g) > tug > 0. Thus the anyon 
energies e A (ki,g) also increase monotonically with i. We therefore have, 

N(e ki ,g)=N(k h 0) = i-l (9) 

To get some feel for the anyon number densities, let us consider the ground state of the N 
particle system. The set {ki} is (1,2,3,...,N). The corresponding set of shifted energies is 
(u, (1 + g)u, (1 + 2g)u, (1 + (N — l)g)uj). The density of particles at energy e is then, 

nA{t,g) = - , e<e F = gp 
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= , e>e F (10) 

Where p = ujN is the average density. This is the same as the zero temperature limit of the 
distribution function in eq. ([!]). 

Next we compute how the anyon density changes with g. From eqs.([|) and (0) we see 
that as g increases the particles move to the right in e space with a 'velocity' given by 
uN(e,g). Thus the number of particles crossing a point e when g increases by Ag, is given 
by the velocity at the point of crossing multiplied by the density at that point. We thus 
obtain the differential equation, 

where p(e, g) = uN(e, g). Denoting the thermal average of p(e, g) as p(e, g) and using the fact 
that in the thermodynamic limit, (p 2 (e,g)) = (p(e,g)) 2 we obtain the differential equation, 



WMl-^WhSl. (12) 

The distribution function n,A{e,g) can be obtained from the solution of eq.(|T2|) by using 
n>A(e,g) = d ^gf ■ Thus we need the solution to eq. (|Ti2|) , with the boundary condition, 



9p{e,g) 



' (13) 



de 

We will now show that eq. ([12|) along with the boundary condition ([13]) is satisfied by , 

p(e,g)= f e de' \ , (14) 
Jo w{e',g) + g 

where w(e,g) is determined through eq.@. By changing variables from e to w in eq.(|l4|) 
the integral can be done to get, 

/ \ 1/, / w(e) . , , w(0) 

^«» = ^ (i "(TT^) ) - ( " ( TT^5) )) (15) 

The condition p(oo,g) = p and the fact that lim e ^ 00 w(e, g) — > oo implies that w(0) is 
independent of g and e _/3/1 = w(0) 9 (l + u;(0)) 1-9 . It can be verified that, 

Tpw e = -7: r- — r P (e, g ), 16 

a 1 + ^ 

^ w(e) = /3^TT^)' (17) 

Using the above results it can be verified that the form of p{e,g) in eq.(|nj) does satisfy 
eq. flT2D . Thus the anyon distribution function in the CSM is given by, 

n A (e,g) = —, — r- (18) 
w(e,g) +g 

which is exactly the distribution function derived by Wu @ for the case of constant spectral 
distribution. The additivity property of the anyon energies suggests that the grand partition 
function should be expressible in a factorized form. If we assume that the usual definition 
of the distribution function for the average number of particles holds for anyons, namely, 

d 

n A (e, g) = —lnz(e) (19) 



where lnZ A = J °° de lnz(e). Then it follows that we must have, 

POD 

lnZ A = / de ln(l + w~\e)) (20) 
Jo 

If this were true, it would show that the particles in CSM have the basic properties of an 
ideal gas. We now show that this is indeed the case. 

The N-particle partition function for the spectrum in eq.(§) is given by 

Zi = JW^ZS, (21) 

where (3 = (3uj and Z^ is the N particle fermion partition function. Setting g = 0, the 
bosonic partition function is obtained, 

Z> = je^zi, (22) 

Combining eqs.fl21|) and (p2[) we may write the anyon partition function as, 

Z A = {Z F N y{Z B N y~*. (23) 
The grandpartition function may also be written in the form, 

oo oo 

Z A = £ eP^Zi = £ (e^(Z£)) V^^*)) 1 ^ (24) 

N=0 N=0 

where we have set the anyon chemical potential = gyLp+il—g)^- in the thermodynamic 

1 dlnZ A 
P 9/J.a 



limit, the sum is saturated at the value of N = N where N = 4 9inZ A g ran( j partition 



function can then be written as, 

Z A (f3, = (Z F (P, ^)y(Z B (P, ^ B )f-\ (25) 

where yLp and /jb are determined by the conditions, N = 4 g = ^ g ■ The solution 
to these conditions is e -^F = w (0),e'^ B = l + w(0). It can be verified that this ensures 
that = N. We therefore have 

P OUA 

roc roc 

In Z A (p, n)=g / dk ln{l + e^'^) - (1 - g) / - S*-"^) (26) 

Jo Jo 



We now make the change of variables (1 + e^ k = 1 + w 1 in the first term and (1 
e -/3(fc-Ms)) = 1 — (1 + w)^ 1 in the second term. We then have, 



oo 



lnZ A ((3,p) = I" dw Z + Q M l + w- 1 ) (27) 

Jw(0) W[l + W) 

We now again make a change of variable from w to e using eq.@. The partition function 
then gets written as, 

lnZ A = j deln(l + w-\e)) (28) 

•J 

This is exactly the factorized form in eq . (p0|) . 

We now consider the equation of state of the ideal gas of anyons in CSM. This has already 
been done by Sutherland [[□]] who however chose to work with the fugacity expansion. The 
coefficients of terms in this expansion in general depend on g and appear quite complicated. 
However, the coefficients of the virial expansion are extremely simple. From the product 
form of the partition function in eq. (p5|), the virial coefficients can all be computed exactly 
since the virial expansion for the free Fermi and Bose gas is exactly computable when the 
density of states is constant. We have, 

oo 

lnZ A = (3P = p J> (Pp) 1 - 1 (29) 
l=i 

where bi are the virial coefficients for the ideal Fermi or Bose gas for I ^ 2 ( bf = bf V I ^ 
2 for constant density of states) and bi = — ^(| — g). It is interesting to note that the 
interaction affects only the second virial coefficient. As we had shown it is the second 
virial coefficient that determines the exclusion statistics. Thus the -4 interaction in this 
system modifies the equation of state in a minimal way. It can therefore be thought of as a 
purely statistical interaction in this system. 

In conclusion, we have shown by applying our definition [|TJ that in the CSM, the particles 
have fractional exclusion statistics. We have then constructed the anyon number densities 
whose average is the ideal anyon distribution function discovered in the context of systems 
with a flat dispersion 0. We have also shown that the partition function can be written 
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in the factorized form characteristic of an ideal gas. Finally we have computed the virial 
expansion and have shown that only the second virial coefficient is affected by the interaction. 
The system thus possesses all the properties we expect of an ideal gas. Thus we have shown 
the system to be an ideal gas with fractional statistics. 
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